Evaporative Cooling of a Guided Rubidium Atomic Beam 
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We report on our recent progress in the manipulation and cooling of a magnetically guided, high 
flux beam of 87 Rb atoms. Typically 7x 10 9 atoms per second propagate in a magnetic guide providing 
a transverse gradient of 800 G/cm, with a temperature ~ 550 /iK, at an initial velocity of 90 cm/s. 
The atoms are subsequently slowed down to ~ 60 cm/s using an upward slope. The relatively high 
collision rate (5 s _1 ) allows us to start forced evaporative cooling of the beam, leading to a reduction 
of the beam temperature by a factor of 4, and a ten-fold increase of the on-axis phase-space density. 

PACS numbers: 32.80.Pj,03.75.Pp 
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I. INTRODUCTION 

Evaporative cooling of atomic clouds confined in mag- 
netic or optical traps has allowed for the realization of 
Bose-Einstein condensation in dilute gases [[J. Coherent 
streams of matter waves have been extracted from such 
condensateSjleading to the achievement of pulsed "atom 
lasers" 0, S, 0, 13 • Nevertheless, the very low duty cycles 
and mean fluxes of the atom laser experiments realized 
to date are crucial inconveniences in view of practical 
applications, such as metrology or nanolithography ex- 
periments. Such applications would indeed require the 
equivalent, for matter waves, of a high flux, continuous- 
wave laser. 

A possible way to realize such a cw atom laser consists 
in outcoupling atoms from a trap periodically replenished 
with new condensates |6(. Another way to achieve this 
goal has been studied theoretically in and transposes 
the technique of evaporative cooling to an atomic beam. 
A high flux, slow and cold atomic beam, transversely con- 
fined by a magnetic guide, can be further cooled by ap- 
plying an energy-selective "knife" which removes atoms 
having a transverse energy above the average. As the 
remaining particles propagate further downstream and 
rethermalize through elastic collisions, the beam temper- 
ature decreases. By lowering the knife as the beam gets 
colder, forced evaporative cooling takes place, and the 
phase-space density of the beam increases. 

The challenging prerequisite in order to apply such 
an evaporative cooling technique is to realize a guided 
atomic beam in the collisional regime, as demonstrated 
for the first time in Here, we report on several im- 
provements of the setup allowing us to achieve a ten-fold 
increase of the beam phase-space density by evaporative 
cooling. An upward slope is implemented along the first 
1.7 m of the guide, leading to a significant decrease of 
the beam velocity. This, in combination with an opti- 
mization of the pulsed injection of atoms into the guide, 
yields a significant increase in the collision rate. We have 
also developed a new detection scheme using absorption 
spectroscopy on an open transition, which improves the 
accuracy of our temperature measurements. 

The article is organized as follows. In section [H] we 
describe the experimental setup, with a special emphasis 



on the implementation of the tilted guide which allows 
us to decrease the beam mean velocity. In section IIIII 
we show the occurrence of non-adiabatic spin-flip losses 
due to the magnetic field vanishing on the guide axis. 
We explain how this can be circumvented by applying 
a longitudinal bias field. We also discuss the necessary 
modification of the analysis made in [i| of the tempera- 
ture measurements deduced from RF spectroscopy on the 
atomic beam. Finally, in the last section, we report on 
what is, to our knowledge, the first significant gain in the 
phase-space density of a guided atomic beam, obtained 
by continuous evaporative cooling. 



II. EXPERIMENTAL SETUP 
A. Magnetic guide 

The magnetic guiding [13, IH IH El Q E3 of our 

atomic beam is ensured by a two-dimensional quadrupole 
guide. The guide axis being chosen as the z-axis, the 
magnetic field reads (bx,—by,§), where 6 is the mag- 
netic field gradient. The field increases linearly with the 
distance from the axis, where it vanishes, which leads 
to non-adiabatic spin- flip losses (see section ^Qj. Con- 
trary to what happens with a spherical quadrupole field, 
this loss mechanism can be circumvented very simply by 
adding a longitudinal bias field B$. The magnetic field is 
then non-zero everywhere. Atoms prepared in a low-field 
seeking state with a magnetic dipole moment /x thus ex- 
perience a semi-linear cylindrically symmetric potential 
given by: 



U(r) = fiJB 2 + b 2 r 



(1) 



where r = (x 2 +y 2 ) 1 ^ 2 is the radial distance from the axis. 
In our experiments, we use 87 Rb atoms polarized in the 
\F = 1, niF = —1) state, which has a magnetic dipole 
moment /i = /l«b/2, where /iB is the Bohr magneton. 

We denote by $ the flux of the atomic beam propagat- 
ing in the guide, v its mean velocity and T its temper- 
ature. Depending on the value of the dimensionless pa- 
rameter a = (iBo/(k-QT), the potential felt by the atoms 
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FIG. 1: Simplified sketch (not to scale) of the magnetic guide. The first 1.7 m of the guide is slightly tilted upwards by 
ho = 22 ± 1 mm, while the remaining part is horizontal within ±1 mm. The guide sits inside glass tubes so that radio-frequency 
(RF) waves can be shone on the atomic beam. Two probe beams, one just after the slope, the other close to the end of the 
guide, can be used to detect the atoms. The beam in the Probing zone 1 can also be used as a marking beam for longitudinal 
time of flight measurements (see text). The ellipses numbered from 1 to 10 represent the RF antennas used for the evaporative 
cooling of the beam; the one labelled T represents the RF antenna used to measure the beam temperature. 



is approximately linear (if a <C 1) or approximately har- 
monic (a 3> 1). For arbitrary a, exact analytical ex- 
pressions for the elastic collision rate 7 within the beam, 
as well as for the on-axis phase space density /?, can be 
derived for a semi-linear potential: 

a_ l + 2a £ / jb_\ 2 /far 
1 2tt 3 /2 (1 + a) 2 v \k B Tj V m ' { ' 

= 1 1 3> / fib \ 2 h 3 
9 2tt 1 + a v \k B T J (27rmfc B T) 3 / 2 ' U 

where m is the atomic mass. In Eq. (0), the elas- 
tic collision cross-section a is assumed to be purely s- 
wave (and thus isotropic) and velocity-independent (a — 
7.6 x 10 -16 m 2 ). Both assumptions are known to be valid 
only if T < 100 /iK for 87 Rb. However the temperatures 
investigated in this article range from 150 /iK to 600 /iK. 
For those temperatures, the e?-wave contribution has to 
be taken into account in the calculation of both the col- 
lision rate and the thermalization rate. The inclusion of 
the d-wave contribution in the estimation of the thermal- 
ization rate is not a straightforward task. One problem 
is the anisotropy of the (i-wave differential cross-section, 
which reduces the efficiency of the thermalization for a 
given number of collisions. Another difficulty lies in in- 
terferences between s and d-waves contributions. For the 
particular case of 87 Rb, the (i-wave cross section exhibits 
a resonance [lGI which increases the collision rate around 
300 /iK. This compensates partially for the reduced effi- 
ciency in rcthermalization of d-wave scattering, resulting 
in a nearly constant thermalization rate |17|. Therefore, 
in the following, we use the expression J2J) as a guideline 
to estimate the number of collisions. 

We now turn to the experimental realization of such 
a guide 0]. Currents up to J = 400 A are running 
in four 4.5 m long parallel copper tubes sitting at the 
vertices of a square of side a — 8 mm. The currents 
change sign between adjacent tubes. A ring-shaped, hol- 
low, non-magnetic stainless-steel piece, brazed on the 



copper tubes, is used at the entrance of the guide in order 
to recombine the currents and the cooling water used to 
dissipate the ~ 5 kW generated by Ohmic heating. The 
copper tubes are held together using ceramic pieces every 
30 cm, and the guide rests on other ceramic pieces which 
themselves rest on the inner side of the vacuum chamber. 
For this configuration of currents, the magnetic held gra- 
dient is b = A/iqI /(ira 2 ), and can therefore reach a max- 
imum value of 1 kG/cm for I = 400 A. The magnetic 
field gradient is actually smaller by a factor of three at 
the entrance of the guide (where a = 14 mm) and reaches 
its asymptotic value only after 40 cm as the distance a is 
progressively reduced. The goal of this tapered section is 
to compress adiabatically the beam during its entrance 
into the guide, in order to maximize the collision rate 7. 

The whole guide sits in an ultra-high vacuum cham- 
ber (see Fig. ^| consisting essentially of two glass tubes 
allowing us to use RF antennas located outside of the 
vacuum chamber, since metallic tubes would screen the 
oscillating magnetic field. Two 25 1/s ion pumps located 
at z — 1.7 m and z = 4.3 m from the guide entrance en- 
sure that the pressure remains around 10~ n Torr at the 
pumps locations; however, due to the poor conductance 
of the glass tubes, the effective lifetime in the second 
part of the guide, measured by monitoring the flux at 
z = 1.7 m and z = 4.3 m for various beam velocities, is 
limited to ~ 15 s. We find experimentally that cooling 
the vacuum chamber to below —50° C increases the flux 
measured at the end of the guide by 25% to 50%. For 
this purpose, we use a flow of cold nitrogen gas emanat- 
ing from a liquid nitrogen dewar. This increase is a clear 
indication that the main source of background gas is the 
outgassing from the glass walls. In addition, we have 
found that the repeated, fast switching-off of the guide 
current leads to a degradation of the vacuum, presum- 
ably by causing micro-leaks to appear at the level of the 
brazing between the copper tubes and the recombination 
piece at the entrance of the guide. For this reason, we 
do not turn off the magnetic field of the guide when we 
measure the atomic beam flux, and we consequently de- 
veloped a flux measurement technique (see § III Cfl which 
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is almost insensitive to the presence of the magnetic field. 



B. Loading of the guide 

In order to produce a high-flux, low velocity, cold 
atomic beam propagating in the guide, we use a pulsed in- 
jection scheme as described in detail in Ref. We recall 
here briefly the various steps involved in the 180 ms se- 
quence used to inject a "packet" of atoms into the guide. 

A slow (~ 20 m/s) atomic beam, with a flux on 
the order of 2 x 10 11 atoms/s, is produced by decel- 
erating, in a Zeeman slower (not shown in Fig. Q, an 
atomic beam effusing from an oven at 140° C. It is 
used to feed a three-dimensional magneto-optical trap 
(MOT). The MOT magnetic field reads (b x x, b y y, b z z) 
with (b x ,b y ,b z ) = (5.4,-6,0.6) G/cm. This anisotropy 
in the held gradients results in an elongated shape of 
the cloud. After 100 ms of loading, the MOT contains 
about 1.8 x 10 9 atoms, and has a transverse rms radius 
of about Ar = 1 mm and a length of 35 mm along z. A 
3 ms launching phase occurs, in which the MOT mag- 
netic fields are switched off, and a small frequency off- 
set 5v is applied between the front and rear pairs of the 
MOT beams lying in the horizontal plane. In this mov- 
ing molasses configuration 19], the cooling occurs in a 
frame moving with a velocity u; n j oc 8 v. This way, the 
atom cloud is set in motion towards the guide entrance, 
with an injection velocity Vj n j adjustable between 30 and 
250 cm/s. In the last 1.5 ms of this phase, the atoms are 
further cooled by increasing the laser detuning from the 
cooling transition to —72 MHz. This results in a final 
cloud temperature around 40 /iK, hence a rms velocity 
Av = 6 cm/s. 

A 700 /is optical pumping stage is used to pump the 
atoms into the \F = 1, m F = —1) ground state, with 
an efficiency of ~ 70%. A two-dimensional quadrupole 
magnetic "pre-guide" is then turned on for 80 ms, with 
a transverse gradient in the x-y plane of 85 G/cm. It 
prevents the atom cloud from expanding transversally 
and from falling due to gravity while it propagates to- 
wards the guide entrance. In order to "mode-match" the 
pre-guide stiffness with the cloud size and temperature, 
one can in principle superimpose a longitudinal bias field 
B\. This makes the pre-guide harmonic, with an angu- 

— 1/2 

lar frequency ui\ o: B 1 . For the parameters of our 
atomic cloud, the field needed to obtain lo\ = Av/Ar 
is B\ = 60 G. Nevertheless, we find experimentally that 
applying such a bias field does not lead to appreciable 
changes in the beam temperature and flux. Moreover, 
the "center matching" (i.e., the overlapping of the atom 
cloud axis with the guide axis z), which is a crucial step 
in order to achieve a low temperature of the guided beam, 
is greatly simplified in the absence of bias field as there 
is no gravitational sag, whereas a bias field B\ — 60 G 
leads to a 3 mm sag. 

The atomic packets containing 1.3 x 10 9 particles, and 
injected with a rate of (180 ms)' 1 , expand longitudinally 



(due to the dispersion of longitudinal velocities) while 
they enter into the guide, and finally overlap, resulting, 
after a propagation over 50 cm, in a continuous beam 
with a flux $ = 7 x 10 9 atoms/s and a mean velocity 

The timings cited here (100 ms of MOT loading and 
80 ms of pre-guiding) are optimized for Vj n j = 90 cm/s. 
For other injection velocities, these values have to be ad- 
justed accordingly. However, it is important to point 
out that the pulsed injection scheme used here prevents 
us from injecting the beam at arbitrary low velocities. 
Indeed, the pre-guiding time r pg must be long enough 
to let the packet propagate over a length I of typically 
10 cm so that the preparation of the next packet does not 
affect the previous one significantly. Therefore, for injec- 
tion velocities much lower than 90 cm/s, the pre-guiding 
time T pg ~ ^/finj would become prohibitively long. This 
would lead to a low repetition rate affecting the value of 
the coupled flux <E>, since the number of atoms that one 
can accommodate in the MOT quickly saturates. 



C. Detection of the atomic beam 

Absorption spectroscopy is used to monitor the beam 
flux, using two laser beams located in the middle and at 
the end of the guide (Probing zones 1 and 2 on Fig.^). 

We have investigated two detection schemes. The first 
one consists in using probe beams resonant with the 
|5 2 S 1/2 ,F = 2) -> |5 2 P 3/2 ,F' = 3) transition, with a 
small amount of repumping light used to first repump the 
atoms from |5 2 Si /2 ,F = l,m F = -1) to |5 2 Si /2 , ^ = 2). 
This has the advantage of yielding a continuous signal if 
the probe beam is locked on the transition. This allows, 
for example, to monitor the time dependence of the flux. 
The velocity measurements described in the next section 
are realized this way. 

However, the guide magnetic field is present during 
those measurements. It leads to a large inhomogeneous 
broadening of the resonance line, largely dependent on 
the field gradient b and on the probe beam polarization. 
We find that in practice, this measurement method is not 
reliable if one wants to measure the atomic beam flux in- 
dependently of the transverse energy of the atoms. For 
the temperature measurements through RF spectroscopy 
(see below), where one is interested in such a measure- 
ment, this method can thus lead to erroneous results. 

We therefore developed another technique, consisting 
in monitoring the absorption of the repumping beam 
alone. The atomic transition is obviously an open one, 
and the absorption is therefore quite small, on the 1% 
level. Nevertheless, provided the light intensity is high 
enough so that any atom, regardless of its position or 
velocity, is addressed by the laser, the number of pho- 
tons absorbed per second is proportional to the atomic 
flux. In particular, for the transition |5 2 S 1 / 2 ,F = 1) — > 
|5 2 P 3 / 2 , F 1 = 2) in which no dark state exists (in contrast 
to what happens for F' = or F' = 1), each atom scat- 
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FIG. 2: Velocity v of the beam after the slope as a func- 
tion of the injection velocity. The magnetic gradient is 

b = 800 G/cm. The solid line represents y v f n j 2gho with 

ho = 22 ± 1 mm (the shaded area represents the uncertainty 
in v due to the uncertainty in ho). 

ters on average two photons before being pumped into 
the \5 2 Si/2,F = 2) level. This is therefore a convenient 
way to measure the absolute atomic flux in the presence 
of the inhomogeneous magnetic field of the guide. We 
find a good agreement between the flux deduced by this 
method and the one inferred from the loading rate of the 
MOT. 

In practice, for measurements of relative flux, we ob- 
tain the better signal to noise ratio when we scan 70 times 
per second the laser frequency across the \5 2 Si/ 2 ,F = 
1) — » |5 2 P 3 / 2 , F 1 = 0, 1, 2) transitions, at a rate of about 
600 MHz/ms. For a beam velocity of v = 60 cm/s, the 
atoms which are in the probing region (of length 1 mm) 
for a given scan were thus located ~ 1 cm upstream dur- 
ing the previous scan and they were not affected by the 
probe. We find that, in contrast to the previous detec- 
tion scheme, the spectra are remarkably insensitive to b 
and to probe polarization. This robustness of the flux 
signal thus leads to a much more accurate temperature 
determination. 



D. Implementation of a slope to slow down the 
beam 

The number of collisions N co \ experienced by an atom 
during its propagation along the whole guide scales as 
1 jv where v is the beam mean velocity. Since the value 
of N co \ determines the potential gain in phase-space den- 
sity that one can achieve through evaporative cooling [9( , 
one needs to minimize v. But as pointed out in § 111 Bl the 
injection velocity Wj n j cannot be arbitraryly low because 
of the pulsed character of our loading process. There- 
fore, exploiting the relative mechanical flexibility of the 
guide and the presence of bellows in the vacuum system 
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FIG. 3: Flux $ measured at the end of the guide as a function 
of the injection velocity Ui n j. The magnetic guide gradient is 
b — 800 G/cm. The timings of the launching sequence are 
kept constant with the values of Section III Bl One clearly 
sees that for low injection velocities, the flux decreases very 
significantly. 

(see Fig. H), we implemented an upwards slope on the 
first 1.7 m of the guide, the height of the guide rising 
by ho = 22 ± 1 mm. The remaining part of the guide is 
horizontal within ±1 mm (see Fig. 0. The purpose of 
such a slope is to decrease the mean velocity of the beam 
and thus increase iV C oi. In the absence of collisions, the 
product of the mean velocity of the beam by the longi- 
tudinal velocity dispersion remains constant in the slope 
as a result of Liouville's theorem. In a collisional beam 
this effect yields a slight increase of temperature by an 
amount 5T/T ~ (2/7)Sv/v [U- 

We have performed velocity measurements in the sec- 
ond part of the guide with a longitudinal time-of-flight 
technique. The laser beam in the Probing zone 1 on 
Fig. ^ located at the position z\, is scanned at a 200 Hz 
rate across the transitions |5 2 Si/2,-F = = —1) — > 

|5 2 P3/2, F', mpi) and is used as a marking beam. For 
a power of 6 /xW, typically 30% of the atoms are re- 
moved from the beam. At time t = 0, we extinguish the 
marking beam, and we monitor with the probe located in 
the Probing zone 2 (locked on the cycling transition) the 
flux arriving at z<i . The time-dependent signal is fitted 
according to 

where erf is the error function, with T and v as adjustable 
parameters. 

Fig. [21 represents the final velocity v measured this 
way, as a function of the injection velocity Vi n j. One 
can see a clear reduction of the velocity, with v close to 

the velocity yj if n j — 2gh expected for a single particle 

having an initial velocity v m j and "climbing" a height ho 
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(solid line). 

The use of a slope to slow down the beam is neverthe- 
less limited by the fact that one must stay in the super- 
sonic regime: the Mach number, defined as M. = v/Av, 
where Av is the longitudinal velocity dispersion, must 
remain significantly larger than 1. If this is not true, 
many particles would move backwards, thereby reduc- 
ing the flux (see Fig. |3J) and increasing the incoming 
particles energy through collisions with a high relative 
momentum. In practice, we find that for v m j = 90 cm/s, 
the final velocity v ~ 60 cm/s (see Fig. [2J is still high 
enough for the beam to stay sufficiently supersonic (for 
T ~ 600 fiK, Av ~ 25 cm/s). The net gain, induced 
by the use of a slope, on the number of collisions N co \ is 
typically a factor of two for this choice of parameters. 
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III. FROM LINEAR TO HARMONIC GUIDING 
POTENTIAL 

A. Spin-flip losses 

In a purely two-dimensional quadrupolc guide, the 
magnetic field vanishes along the guide axis z. There- 
fore, atoms having a low angular momentum experience 
a rapid variation of the direction of the magnetic field 
when they approach the axis, and can undergo a non- 
adiabatic spnVfhp transition to an untrapped state & 
This loss mechanism is more pronounced if the beam en- 
ters the collisional regime, since collisions redistribute an- 
gular momentum between atoms and continuously lead 
to new losses. Adding a longitudinal bias field to the two- 
dimensional quadrupolc leads to a semi-linear transverse 
trapping potential JQ, in which the magnetic field never 
vanishes, which can strongly suppress those Majorana 
losses. 

To study the occurrence of such spin-flip losses in our 
beam, we use a solenoid wound around the guide, with 
a diameter varying from 30 to 50 cm and 20 turns per 
meter, yielding a longitudinal bias field Bq = 0.25 G/A 
along z. The resulting field is uniform within ±10% along 
the whole guide length. 

When we apply a sufficiently large bias field (B > 1 G) 
the flux in the guide is enhanced by as much as 30% 
(see Fig. with respect to its value for Bq = 0. 
The large width of the loss curve (1 G FWHM) indi- 
cates that the observed losses are not, strictly speak- 
ing, due to Majorana spin-flips, but are probably in- 
duced by high-frequency oscillating fields (due to the 
noisy environment) having an amplitude of a few tenths 

of gauss in mm. 

In practice, we operate the magnetic guide with a Bq = 
1 G bias field. This changes slightly the shape of the 
transverse trapping potential. If the beam temperature 
T is large {k^T 3> (J-Bq) the trap is still approximately 
linear; but for lower temperatures the relevant part of 
the potential is harmonic. 



FIG. 4: Measured flux at the end of the guide (normalized to 
the flux in the absence of an applied bias field) as a function 
of the applied bias field ^applied- The curve is not centered on 
because of the presence of a -B oart h ~ 120 mG component 
of the Earth magnetic field along the guide. The bias field is 

Bo = B a ppli c( j — B calt h- 



B. Temperature measurements in a semi-linear 
potential 

The exact shape of the potential needs to be taken into 
account in the temperature measurements that we per- 
form with RF spectroscopy For a linear potential, 
we extract the beam temperature T from the fraction ip 
of atoms remaining in the beam after evaporation by a 
RF antenna operating at frequency v x f. If this frequency 
is small compared to k^T/h, then very few atoms will 
fulfill the RF resonance condition at some point in their 
trajectory and will be evaporated. In a similar way, if 
hi/rf 3> fceT, the resonance condition is fulfilled on a 
radius large compared to the thermal radial size of the 
beam, and almost no atoms are evaporated. In between, 
the function ip goes through a minimum. The same fea- 
ture remains qualitatively true for a semi-linear potential, 
but the precise lineshape of this RF spectrum depends on 
the value of the bias field. 

We therefore developed, for an arbitrary bias field 
Bo, semi-analytical calculations allowing us to extract 
the temperature from the curve depicting the fraction 
y(frf) °f remaining atoms after the temperature mea- 
surement antenna, as a function of its RF frequency 
f r f. The method follows the one used in |9j. The curve 
(p(vif) now depends on two dimcnsionless parameters, 
rj = h{y T i — co)/(^bT) (where vq is the frequency cor- 
responding to the "trap bottom" , satisfying hv§ = (J.Bq) 
and a = (xBq/QcbT). In the limit a — > one has a linear 
trap, leading to the approximate formula 

ip( V ,0) ~ 1 - 1.65 77 113 e- a92 \ (5) 

while for a — > oo one recovers the exact formula for a 
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harmonic confinement: 

(f(n,oo) = 1 - ^/wfje' 71 . (6) 

From our calculation, we find a good interpolation be- 
tween those two limits is provided by the following em- 
pirical formula: 

Vap P rox(r?,a) 1 - 1.7 e-°-9V- 1_0 - 4arctan(3 - 6a) ■ (?) 

We checked the validity of the approximate formula Q 
by generating, with a Monte-Carlo simulation, ensem- 
bles of atoms at thermal equilibrium in the semi-linear 
potential, simulating the curve (p(rj, a) by applying the 
evaporation criterion, and fitting this curve by Eq. 
We find that the temperature obtained in this way is reli- 
able within ±5% for a range of variation of 0.1 ^ a ^ 10. 
The precision with which we can infer temperature ratios 
is actually better. We have checked that, for a given B , 
when two temperatures Tf lm and T| lm are used to obtain 
simulated RF spectra, the ratio of the temperatures T-f * 
and * obtained by fitting those spectra with Eq. (J7J 
differs by typically less than one percent from the actual 
ratio Tf m /Tf im . 

With this fitting function, we investigate the effect of 
spin-flip losses on the beam temperature. In the absence 
of an applied longitudinal field, we observe that the beam 
temperature is ~ 8% higher than when one applies a bias 
field large enough to avoid spin-flip losses. This can be 
easily understood since the atoms lost by spin-flips have 
a low angular momentum, and on average a small energy. 
Therefore their loss leads to a higher beam temperature 
after rethermalization. 



IV. EVAPORATIVE COOLING OF THE BEAM 

It is important to emphasize that applying evaporative 
cooling on an atomic beam is significantly more difficult 
than on a trapped cloud of atoms: (i) the pulsed injec- 
tion results in a longitudinal dilution of the atom packets, 
and the resulting decrease in the atomic density yields a 
low initial collision rate 7, (ii) the time allowed for evap- 
oration is limited by the guide length and by v, and not 
only by the quality of the vacuum, and (iii) the collision 
rate and phase-space density scale less favorably in view 
of evaporative cooling for a two-dimensional confinement 
than for a three-dimensional trap 9]. For all those rea- 
sons, a significant gain in the phase-space density of a 
magnetically guided beam had never been achieved to 
date. 

In this section, after detailing the two evaporation 
methods that we have developed, namely radio- frequency 
and microwave evaporation, we report on the evapora- 
tive cooling of the atomic beam using several evaporation 
zones in order to significantly increase the phase-space 
density. 



A. Radio- frequency evaporation 

Evaporative cooling relying on a position-dependent 
transition to an untrapped state can be easily accom- 
plished using a two-photon radio- frequency transition be- 
tween the Zeeman sublevels \F — l,mp = —1) and 
\F = l,mp = 1). Atoms crossing a cylinder of axis 
z, with a radius i? C vap fulfilling U(R cvap ) = hv T t, where h 
is Planck's constant and u T { the RF frequency, are evap- 
orated. If this radius is higher than the beam's thermal 
size, essentially energetic atoms are removed, and after 
rethermalization, the beam gets colder. For temperatures 
in the range of a few hundreds fiK, the corresponding 
frequencies are in the tens of MHz. The RF antennas 
that we use consist of a single loop with a 50 resistor 
in series, in which we couple a power of 27 dBm. This 
corresponds to an amplitude of the oscillating magnetic 
field, at the level of the atomic beam, of 50 mG. 

We can measure the spatial range over which an an- 
tenna is efficient by switching on a RF antenna for a 
short time (a few tenths of milliseconds) and monitor- 
ing the temporal width of the dip appearing in the flux 
signal. We find that the typical range is around 20 cm. 

Nevertheless, for some specific frequencies (especially 
around 40 MHz) we observe that the efficiency range 
of the antenna considerably increases, essentially reach- 
ing the whole guide length. Those resonance frequen- 
cies depend on the position z ant of the antenna on the 
guide. This fact, added to the similarity between the 
guide length and the wavelengths corresponding to the 
resonance frequencies, leads us to suppose that the RF 
antennas can couple to the guide, in which they induce 
currents, making it act as an antenna itself. 

Although those resonances are well localized in the fre- 
quency domain, and therefore do not prevent the use of 
RF evaporation for the purpose of evaporative cooling, 
we tried to circumvent this problem with evaporation in 
the microwave domain by inducing transitions between 
different hyperfine states. 

B. Microwave evaporation 

An electromagnetic wave resonant with the hyper- 
fine transition |5 2 S 1/2 , F = 1) -> |5 2 S 1/2 , F = 2) at 
i^hf = 6834.7 GHz can also be used to selectively remove 
atoms from the beam 25]. Indeed, a magnetic dipole 
transition from the trapped state \F = l,mj? = — 1) 
can only excite atoms to the anti-trapped states \F — 
2,m,F = — 2), \F = 2, 771 p = —1), or to the untrapped 
state \F = 2,777,77 = 0). Contrary to the RF-case, this 
gives three distinct evaporation radii. This leads actually 
to a better efficiency of the evaporative cooling process 
compared to single-radius evaporation, as was discussed 
in |9j for the case of evaporation with a continuum of 
radii spanning the range [i? cva p,oo). 

We calculate the fraction ^(!/ mw ) of atoms remaining 
after such a microwave antenna operating at z/ mw , in the 
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FIG. 5: (a): Fraction of remaining atoms after the microwave 
horn antenna as a function of 77. The experimental points are 
in good agreement with theory (solid line) ; the points at 77 — 6 
and 8 are slightly above the expected value, probably because 
of a lack of microwave power to reach a 100% efficiency in the 
evaporation, (b): Temperature change induced in the beam by 
evaporation with a microwave at the parameter 77, followed by 
rethermalization. The solid line is the theoretical curve; the 
points are experimental results of evaporation performed in 
the interesting region where the phase-space density increases 
significantly. For 77 — 8, the phase-space density of the beam 
increases by 2.6. For those experiments, the mean velocity is 
v = 0.8 m/s and the guide gradient 600 G/cm. 

same way as for the single-radius evaporation scheme. 
Defining 77 as h(vh{ — I'mw) /(ksT), we work out the fol- 
lowing approximate formula, valid for a linear potential: 
ip(ri) = 1 - 0.89 77°- 67 e-°- 27 ''. The shape of this curve 
plotted on Fig. (a) is qualitatively similar to the one 
obtained with RF evaporation, but with a much larger 
width and a more pronounced minimum. 

In practice, microwave evaporation is implemented by 
using a 10 dB gain horn antenna powered by a microwave 
synthesizer delivering typically 20 dBm after amplifica- 
tion. The directivity of such an antenna allows us to 
use the microwave evaporation scheme even in locations 
where nearby metallic walls prevent the use of RF evap- 
oration. For instance we can use it in the metal cross 
located in the middle of the guide, by shining the mi- 
crowave through optical viewports (see Fig. ^| . We have 
used this evaporation scheme in order to measure the 
beam temperature and found a very good agreement (less 
than 10% difference) between microwave and RF spectro- 
scopic measurements of temperature. 

We have also calculated the temperature change T' /T 
experienced by the beam upon evaporation with a mi- 
crowave antenna followed by rethermalization. Fig. [31(b) 
shows the corresponding theoretical curve, and the ex- 
perimental data points obtained for evaporation with 
77 = 6, 8, 10 and 12. The beam temperature was mea- 
sured with a RF antenna. The agreement is excellent, the 
theoretical curve having no adjustable parameters. For 



the point corresponding to 77 = 8, the gain in phase-space 
density is 2.6, which is significantly higher than the max- 
imum gain of 1.9 that one can obtain with (single-radius) 
RF evaporation. 



C. Ten-fold increase of the beam phase-space 
density with several antennas 

In order to achieve a significant gain in phase-space 
density, we use the following setup (Fig. We choose 
«inj = 90 cm/s and b = 800 G/cm. Five RF antennas, 
separated by ~ 20 cm, are placed between 0.5 ^ z ^ 
1.3 m with frequencies decreasing from 49 to 25 MHz, 
then the horn antenna at z — 1.7 m is operated at 
c mw = 6774 MHz, and finally five more RF antennas 
are placed at 2 z ^ 2.8 m with frequencies decreasing 
from 21 to 14 MHz. We use such a scheme, in which ob- 
viously no complete rethermalization can occur between 
successive RF antennas 1 , because it is less sensitive to 
a potential lack of efficiency of an antenna for some fre- 
quency. In addition, this scheme is closer to the usual 
evaporation ramp in the time domain used for evapora- 
tive cooling of trapped atoms than the step-like evap- 
oration process that would consist in letting the beam 
rethermalize completely after each antenna. 

The last evaporation antenna being placed at z = 
2.8 m, and the analysis antenna used for the tempera- 
ture measurements at z — 4 m, the beam can rether- 
malize completely after the last antenna, and therefore 
we measure the temperature of a beam which is really in 
thermal equilibrium. 

Fig. 0shows the RF spectra from which we deduce the 
initial (without evaporation) and final (after evaporation 
by the 10 RF and microwave antennas) temperatures 1] 
and T f . We find 7] = 574 ± 9 (iK, and T f = 164 ± 6 (jK. 
The corresponding flux reduction is $f/$i = 0.13 ±0.02. 
From this we infer a gain in phase-space density by a 
factor 10.4j/J;J. The spectrum giving Tf corresponds to 
a well thermalized beam (the spectrum does not reach 
the value 1 for, e.g., v T { — 14 MHz, which would be 
the signature @ of the absence of rethermalization after 
the last evaporation antenna). This configuration cor- 
responds to the existence of a steady-state temperature 
gradient dT/dz ~ 200 /iK/m along the beam. 

Those measurements require a good stability of the 
incident atomic flux. Indeed, a small relative variation 
<5<&i / <E»i can be amplified very significantly since a de- 
crease in $; will result in less collisions between succes- 
sive antennas, and thus in a more pronounced reduction 
of the flux than in the collisional case. We observed that 
fluctuations in $; by 30% could lead to fluctuations above 



1 Even at the low velocity achieved in the second, horizontal part 
of the guide, on average only ~ 1.5 collisions per atom can occur 
over 20 cm. 
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FIG. 6: Radio-frequency spectra used to determine the tem- 
perature of the atomic beam with and without evaporation by 
10 RF antennas and the horn antenna. The temperature re- 
duction from 574 to 164 /iK, accompanied by a flux reduction 
by a factor 7.7, leads to a ten-fold increase of the phase-space 
density. The arrows represent the evaporation frequencies of 
the 10 RF antennas. The two longer arrows correspond to 
RF frequencies that would lead to the existence of evapora- 
tion radii coinciding with the ones resulting from microwave 
evaporation by the horn; another one at 60.7 MHz is not 
shown on the figure. 

50% in $f. For the data of Fig. we therefore checked 
that the absolute flux was almost constant and equal to 
its maximal value during the whole data acquisition. 



V. CONCLUSION AND PROSPECTS 

In this article, we have demonstrated the first signif- 
icant gain in the phase-space density of a magnetically 
guided atomic beam by means of continuous evaporative 
cooling. This ten-fold increase brings the beam's phase- 
space density to about 2 x 1CP 7 . In order to gain the 
seven remaining orders of magnitude required to reach 
quantum degeneracy, the temperature needs to be de- 
creased so much that the guide confinement will become 
harmonic. As pointed out in [jj, the evaporation becomes 
less efficient for this confinement as compared to the lin- 
ear case, and the elastic collision rate cannot increase 
significantly during the evaporation process. 

In order to have a quantitative picture of the whole 
evaporation ramp all the way to quantum degeneracy, we 
therefore calculated for arbitrary values of a, the change 
in flux and temperature for a given value of rj. This 
allows one to calculate, through Eq. (J2J and ©, the gain 
in phase-space density and the evolution of the collision 
rate. 

The resulting "evaporation trajectories" followed by 
the beam in the plane (T, when it crosses successive 
single-radius evaporation antennas are plotted on Fig. [7| 




Temperature T (/iK) 

FIG. 7: Theoretical "trajectories" to quantum degeneracy fol- 
lowed by the beam in the plane (T, $) during evaporation 
with many antennas, at constant rj. The grey scale represents 
the on-axis phase-space density p (which scales as $/T 7,/2 in 
the linear regime a <C 1, and as $/T 5,/2 in the harmonic 
regime). The atomic beam propagates with a velocity of 
60 cm/s in a magnetic guide having a transverse gradient of 
800 G/cm, with a longitudinal bias field of 1 G. The initial flux 
is 7 x 10 9 atoms/s, and the initial temperature 570 pK. Each 
symbol in the various evaporation trajectories corresponds to 
the action of one antenna. For higher evaporation parameters 
?], the final flux and temperatures (when p = 1) are higher, 
but this requires more antennas (the number of antennas re- 
quired to reach p = 1 for r\ — 3, 4, 5, 6 and 7, is 41, 64, 106, 
181, and 306, respectively). The black squares represent the 
experimental beam parameters before and after evaporation 
with the method described in § IIV Gl 



and |H1 The frequency of each antenna is adapted so as 
to evaporate at constant r\ parameter, and we assume 
that full rethermalization occurs after each antenna. One 
sees that evaporating at relatively low r\ minimizes the 
number of antennas needed to reach degeneracy, but im- 
plies to have very low final flux and temperatures, in the 
30 nK range. Simultaneously, the collision rate for such 
low evaporation parameters is decreasing significantly in 
the final evaporation stages when the temperature is low 
and the potential harmonic (see Fig. [SJ . On the contrary, 
evaporating at too large r\ requires too many antennas 
and thus a very long guide. 

A promising strategy probably consists in operating 
the antennas with r\ — 4 or 5, requiring a number of 
antennas between 60 and 100. Nonetheless, due to the 
finite range of the antennas, with such a large number 
of evaporation zones one probably needs to abandon the 
simple model of discrete evaporation zones followed by 
rethermalization |2ti[ , and resort to a continuous evapora- 
tion model as the one used in standard BEC experiments 
and considered in Q. The atomic energy distribution is 
then at any time close to a truncated Boltzmann distri- 
bution |27| , whereas we considered here that enough time 
was allowed between two successive antennas to recover 
an exact equilibrium. Another possibility would consist 
in using surface evaporation [28j by approaching solid 
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surfaces close to the beam, on which the most energetic 
atoms would stick: this would allow for the use of many 
very localized evaporation zones. 



Temperature T (jXK) 



FIG. 8: Same as Fig |7| but with a grey scale representing 
the elastic collision rate 7, assuming that the collision cross 
section is purely s-wave, and constant in magnitude (both 
assumptions fail above typically 100 /iK). The cross-over from 
a linear trap (in which 7 ~ <J>/T 3//2 ) to a harmonic trap (where 
7 ~ $/T 1//2 ) is clearly visible when feT becomes on the order 
of fiBo (here, at T ~ 30 ^K since Bo = 1 G). Runaway 
evaporation, ie constant increase of 7, can only occur for large 
values of n; for n = 4, which requires a moderate number 
of antennas, 7 starts by increasing at the beginning of the 
evaporation, when the trap is still linear, but then decreases 
as the guide becomes more and more harmonic. 
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